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Abstract 

A one-way wave equation is an evolution equation in one of the space directions that de- 
scribes (approximately) a wave field. The exact wave field is approximated in a high fre- 
quency, microlocal sense. Here we derive the pseudodifferential one-way wave equation 
for an inhomogeneous acoustic medium using a known factorization argument. We give 
explicitly the two highest order terms, that are necessary for approximating the solution. 
A wave front (singularity) whose propagation velocity has non-zero component in the spe- 
cial direction is correctly described. The equation can't describe singularities propagating 
along turning rays, i.e. rays along which the velocity component in the special direction 
changes sign. We show that incorrectly propagated singularities are suppressed if a suitable 
dissipative term is added to the equation. 
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1 Introduction 



In this paper we consider one-way wave equations for inhomogeneous acoustic 
media in n dimensions, where n > 2. We assume there is a special space direction, 
which we call depth or the vertical direction, with coordinate denoted by z, the other 
directions are called lateral or horizontal and are denoted by x. The time coordinate 
is denoted by t. The medium is described by its slowness (inverse velocity) v = 
u(z,x) and its mass density p = p(z,x). Let U = U(z,x,t) denote the acoustic 
wave field, and F = F(z, x, t) a volume source, then the acoustic equation is given 
by 

n-l 

-p- x v 2 d] + d x1 P~ l d Xi + d zP - l d z \U = F. (1) 

3=1 
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A one-way wave equation is an equation describing only downward propagating 
waves, whose propagation velocity has positive vertical component, or only upward 
propagating waves, with negative vertical component of propagation velocity. As 
we will discuss here, by solving a one-way equation by progressing in depth in 
one direction an approximation to (part of the) wave field is obtained. One-way 
equations are used in applications in geophysics (see e.g. [1]) and ocean acoustics, 
for which many different numerical methods have been developed. References and 
a discussion of some of these can be found in e.g. [4,2]. 

In an inhomogeneous medium reflections occur for waves with wave length com- 
parable to the scale of the medium variations. Only the high-frequency part can be 
expected to be computed using a method progressing in one direction only. How- 
ever, this is not necessarily a disadvantage. In seismic imaging and migration an 
asymptotic limit is implicit, and the absence of reflections is often an advantage. 

The high-frequency (or singular) part of solutions to wave equations is very well un- 
derstood. High frequency waves propagate along rays, curved trajectories in space. 
We use the theory of microlocal analysis about this (see e.g. the books [3,6,9]). 
Thus we consider approximation of the wave field modulo an error that is C°°, 
hence in Fourier space goes to zero when the frequency goes to infinity, faster than 
any negative power of the frequency. The singularities of the function, that is the 
part that does not go to zero rapidly when the frequency becomes large, can be 
localized in position and direction using the wave front set of Hormander (see the 
mentioned references or [5]). For propagating singularities the position and direc- 
tion determine the ray. Using these ideas we describe precisely how the wave field 
is approximated by solving a one-way wave equation. 

Singularities propagating with velocity that is not horizontal can be described by 
a pseudodifferential evolution equation in z. Such equations are obtained from a 
factorization or decoupling argument, see e.g. [8] or the similar treatment in [9, 
section 9.1], or [6, section 23.2]. However, such equations are not defined at the 
point where the ray is tangent to horizontal. Moreover, it is clear that waves prop- 
agated along turning rays cannot be computed by progressing in one direction in 
depth only. In practice it is desirable that such wave fronts are suppressed. 

In this paper we give a pseudodifferential evolution equation in z that includes 
such a suppression. Thus we obtain a model for numerical one-way wave equation 
methods. We give sufficient detail so that the equation can in principle be used as a 
starting point for discretization. 

Our one-way wave equation is obtained in two steps. First we extend the factoriza- 
tion argument of Taylor [8]. We compute the pseudodifferential equation referred to 
above to the two highest orders that are at least needed for a highest order approx- 
imation of the solution (for such a pseudodifferential equation, as for the ordinary 
wave equation, the leading order term describes the behavior of the rays, while the 
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next term describes the amplitudes to highest order). We also allow for a normaliza- 
tion in the definition of the down- and upward propagating parts of the wave field. 
With a suitably chosen normalization the equation becomes unitary, microlocally. 
We then modify this equation in order to suppress any wave fronts propagating 
along turning rays. We show that the solutions to the one-way equation approxi- 
mate microlocally the real solution. 

The precise formulation of these results is the subject of the next section. Sections 3 
and 4 contain the proofs of two theorems. 



2 Pseudodifferential one-way wave equation and approximation of the solu- 
tions 

We first introduce some notation. The Fourier variables corresponding to x, z and t 
will be denoted by £, ( and r, with Fourier transform defined by 

m, t) = J J e-'K-*-^) f(x, t) dx dt. (2) 

A function ip = ip(x, t, £, r) in C°°(]R Tl x R n ) is a symbol of order m, if there are 
constants C a ^ such that 

\d:X,M^ *, & r) I < C Q>/3 (1 + || (£, r) \\) m -W (3) 

(see a text on pseudodifferential operators). Here a, (3 are a multi-indices a = 
(ai, . . . , a„), cxj E Z> , |a| = «! + ... + a n . Associated with such a symbol ip is 
a pseudodifferential operator, that will be denoted by ip(x, t, D X) D t ), and is given 
by 

i>(x,t,D x ,D t )f = (2tt)-"/ iKx,t,Z,T)f(Z,T)j& x+T Q d^dr. (4) 

The set of symbols is denoted by S m (M. n x W a ). We will typically encounter op- 
erators acting in (x, t), depending on z, but independent of t (of convolution type 
in t), i.e. with symbols ip(z, x, r). A symbol is in S~°° if for any N there is a 
constant C such that \d£ t d£ T iJ)\ < C(l + ||(^, t)||)~ w The corresponding operator 
then maps distributions to C°°. We will write ipi ~ ^ 2 if ^ - ^ 2 £ S*" 00 . 

We recall that, if f(y), y E M™ is a distribution, Hormander's wave front set WF(/) 
contains points (y, rf), contained in the cotangent space or phase space W 1 x M n , that 
give positions and directions associated with the singularities of /, see [3, section 
1.3], or [5, section 8.1]. This set is conic, i.e. if (y : rj) E WF(/) then (y : \rj) E 
WF(/) for all A > 0. If T is a conic subset of R n x R n we say that f = g 
microlocally on T if WF(/ - g) n T = 0. 

Denote by p(z, x, (, ^, r) = x)~ l v{z, x) 2 t 2 — p _1 (^ 2 +C 2 ) the principal symbol 
of P. If f7 satisfies PU = 0, then the singularities of U are in the characteristic set 
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given by 

P (z,x,{,S,t)=0. (5) 

They propagate along null bicharacteri sties, curves in the cotangent space contained 
in the set given by (5) that are solutions to the Hamilton vector field of p (see 
[9, Theorem 6.2.1] or [6, section 23.1]). Parameterizing by time, the differential 
equations for the null bicharacteri sties are 

^ = - *(z, x)- 2 r-% ^ = - v(z, x)-\-X (6) 

— = -TV' 1 — ^L=_ TU ^— (7) 

dt dx ' dt dz' 

So singularities with — £ > 0, satisfy > (downgoing) and waves with — ^ < 
satisfy < (upgoing). 

If a point (z, x, £, r) is given with ||£|| < v(z, x)\t\, then there are two solutions ( 
to (5). These will be denote by ±6 = ±b(z, x, £, r), where 

b = -TUy/l - T- 2 U- 2 i 2 . (8) 

The sign is such that ±6 corresponds to rays with ±|| > 0. We define a set of 
points (z, x, £, r) associated with propagation angles with the vertical < 9, by 

I' e = {(z,x,C,r)\r^ and \\u(z, x)- 1 ^ 1 ^ < sm(9)}. (9) 



To obtain the microlocal one-way wave equations we follow the factorization ar- 
gument of Taylor [8]. In this factorization it is assumed that singularities of F and 
therefore also those of U are not in the set given by £ = r = 0, ( ^ 0. By (5) 
there is a constant C (assuming that v is bounded) such that for the propagating 
singularities we have 

lCl<C|r|. (10) 

The decoupling will be done microlocally where the propagation angle 9 is smaller 
than some fixed angle 9 2 . We define a subset of phase space IR n+1 x W l+1 associated 
to such angles by 

Ie 2 = {(z,x,t,C,C,r) | (z,x,£,t) G I' $2 , |C| < C\t\}. (11) 
Let I±fl 2 denote the subsets of I e , 2 with =Fr _1 £ > 0. 

The down- and upgoing components (u+,u-) are obtained from {U,p^ 1 ^) by a 
2x2 pseudodifferential matrix Q = Q(z, x, D x , D t ). We assume that Q is elliptic, 
i.e. there is a matrix pseudodifferential operator W such that WQ ~ QW ~ Id. 
We take the liberty to denote Q^ 1 = W, even though this is not quite correct. We 



4 



also define sources (/+,/_). We have 

:Hi'f)' (;:H-(> - 

We have the following result about one-way wave equations for (u + , u ). The proof 
is the subject of section 3. 

Theorem 1 For suitably chosen Q and B± = B±(z, x, D x , D t ) the equation 

PU = F microlocally on I 02 . (13) 

holds if and only if 

P 0)+ u + = f (d z — iB + (z, x, D x , D t ))u + = / + , microlocally on Ig 2 , and (14) 
Po,- M - = f (dz — iB_(z, x, D x , D t ))u_ = /_, microlocally on I$ 2 . (15) 

The operator B± can be chosen selfadjoint, with B± and Q satisfying 

& +l irl E^-^-j +order(-l), on l' 02 , (16) 

Q(z,x,Z,t)=\ _ t 1 + order , on 1^ 

\isgn(r)p 2 fl 4 -isgn(r)p id* J \~2 ~2j 

(17) 



where a is defined by a(z, x, £, r) = rr) 2 r 2 — ||£|| 2 . //we choose Q such that 
U = u + + that is Q 1} \ — Qi ^ — X then B± satisfies 

n t r \ I, 1 .,-1 x 1 - 1 \ ^.^a ,.dp _-, 

B±(z, x, {, r) = ± [b + ^ g ^ j + - 

+ order(-l), on (18) 



The highest order term in equation (14) determines the Hamilton flow of the singu- 
larities, hence the rays. The zeroeth order term determines the amplitude. Thus both 
terms need to be incorporated for an accurate highest order approximation of the 
solutions using (14). There are two kinds of zeroeth order terms for B±. First the 
term ^ib' 1 YTjZl |r J^- When computed explicitly it is equal to ±u/f • %J' x {y 2 - 
r -2£2)-3/2 This term mates me operator b + \ib~ l Y%Z\ || £h selfadjoint up to 

zeroeth order. Second there is the term |i||a _1 — |iffp _1 in (18), due to the dif- 
ferent normalization of Q. 

Equation (14) is only microlocal. Outside F d2 the symbol B± is not prescribed, but 
we choose it with real principal symbol and smooth (i.e. without the singularity of 
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the square root). A wave front propagating on some turning ray will be propagated 
incorrectly with this equation. To suppress such singularities we introduce a damp- 
ing term given by a pseudodifferential operator C = C(z, x, D X) D t ). The complete 
one way wave equation will be of the form 

P ±U± = (d z - iB±{z, x, D x , D t ) + C(z, x, D x , D t ))u ± = 0. (19) 

It is assumed here that a solution is sought for z > z (for z < z the sign in front of 
C(z, x, D x , D t ) must be changed). We assume the dissipative term is for waves 
propagating with angle smaller than some given angle 9\, 9 X < 9 2 . The operator 
C is also a pseudodifferential operator with homogeneous, non-negative principal 
symbol c(z, x, £, r). We let the order be 1, so that the length scale associated with 
the decay is proportional to wave length, but this is not essential. Its main property 
will be 

c(z,x,£,t) = for(z,x,£,T)eI' 01 (20) 

c{z,x,£,t) > rj(Z 2 + T 2 )? for (z,x,£,t) outside 4, (21) 

where rj is some positive constant. In addition there is the condition that when 
C — 0, then also a number of its derivatives are zero, see the precise formulation 
below. 

We consider the approximation of solutions to the equation (1) with right hand side 
0, given that the solution U has only singularities propagating in the + direction (or 
only in the — direction). So suppose that U satisfies 

PU = for z>z , (22) 

and assume it has only singularities propagating in the + direction at z = z , in 
other words 

WF(U) n{z = zq, r _1 C > 0} = 0. (23) 

We also assume that the singularities of PU at z = z satisfy (10) at z = z , which 
implies that the restriction U\ ZQ is well defined. Let Q + = Q- = Qi,2- The 
approximate solution is then given by Q+u + , where u + is the solution of 



u+ 



Q+^or^L, (24) 



P+u + = 0, for z > z , (25) 



With a point (z , x, t, r) e V e and a time t there are two associated null bichar- 
acteristics corresponding to the two possible values of ( = ±b(z ,x,£,r). They 
can be parameterized by ^-coordinate of the ray as long as the angle of the velocity 
vector is smaller or equal than 9. We let [Z min (z , x, £, t, r, 9), Z ma _ x (zo, x, £, t, r, 9)\ 
be the maximal interval where this is the case, and we denote the bicharacteristic 
by 7±(z, zq, x, t, ^, t). Let c{z, x, ^, r) be the principal symbol of C. It was shown 
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in [7] that in a solution operator for (25) a pseudodifferential factor occurs with 
symbol 

exp ^ - £ c(-y±(z', z, x, t, £, t)) dz'^J . (26) 

This exponential is equal to 1 when the bicharacteristic 7± stays in I 8l , while it is 
exponentially decaying to for (£, r) to infinity if a finite segment between z and 
z is outside the region where c = 0. Therefore we define a subset of W n+1 x M n+1 
that can be reached from depth z , while staying in I e , by 

J±(zo, 0) = {(*, x, t, C, £,r) | =F r _1 C > and Z mhl (z, x, £, r, 6) < z }. (27) 

We will show U can be approximated by Q+u + in the following way 

WF(Q + u + ) c WF(U), (28) 
Q + u + = U on J+ (20, (29) 
= outside J+(^o, ^2)- (30) 

The same is true with u + , Q + and J + replaced by w_, Q_ and J_. 

In [7] additional assumptions on C were made. We first give an example. Define a 
scalar function h(y) that smoothly goes from constant equal to zero at y < 0, to 
being positive at y > by the formula 
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y<0, 
exp(-l/y)/(exp(-l/y) + exp(-l/(l - y))) < y < 1, (31) 

1 y> 1 



Now define for example 

C = w{z,x,i,r)h{u~ l \\r- l i\\ - sin(^)) (32) 

where w is homogeneous of order 1 in (£, r) and bounded below by some constant 
times V'?" 2 + £ 2 - 

In general we make the following assumptions. We assume that C is given by a 
sum C = c + where C* 1 ) is of order 0. We will also write C^°\z, x, t) = 
c(z, x, £, r). We assume that there is an integer L > 2 such that the derivatives of 
order up to L of and of order L — 2 of C^ 1 ^ satisfy the following bounds 

\did:dl r C^(z,x,£,T)\ < C(l + ||(e,r)||)-l^ i ^(l +C (,,.,e,r)) 1 -^. 

j + 2fc + |a| + |/3| <L, (33) 

for some constant C. It was shown in [7] that (32) satisfies this property for any L. 
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The initial value problems for the operators Po,± an d P± that were defined in (14) 
and (19) have well defined solution operators, that we will denote by E 0j ±(z, z ) 
and E±(z, z ). In [7] it was shown that these are related by a pseudodifferential 
operator K± = K±(z, z , x, D x , D t ) with principal symbol (26), such that E± ~ 
K±Eo,± (with K in a class of symbols more general than that given by (3)). With 
the assumptions on C it follows that K has the property 

K±(z, Zq, x, £, t) ~ 1 on J±(z , 0i), (34) 
K±(z, zq, x, £, t) E S°° outside J±(z , 9 2 ), if z — z > 5, (35) 

if 5 > is some small constant. This results in the following theorem, that gives 
sufficient conditions for the approximation property of equations (28), (29) and (30) 
to hold. The proof is given in section 4. 

Theorem 2 Let U — U(z, x, t) satisfy (22) and (23). Let u + be the solution to (24) 
and (25), where Q + , B + are as in Theorem 1, and C satisfies (20) and (33). Then 
there is K, depending on Q + , B + and C, satisfying (34) and (35) such that for 

z > z 

Q + u + = KU + r, (36) 

withr e C°°{]z , oo[xR n ). 
3 Proof of Theorem 1 

The computation of Q and B± is done by writing (1) as a system of first order in z, 
and then transforming this system. Let V be defined by V = P' 1 ^, and A, A p by 

A p = -p-V^x^ + ^^.p- 1 ^., (37) 

j 

A = A 1 = -v(z,x) 2 d 2 t +Y,d 2 Xj . (38) 

3 

The principal symbols of these operators are a(z,x,^,r) = u(z,x) 2 t 2 — ||£|| 2 , 
a p (z,x,^,r) = p^ 1 u(z,x) 2 T 2 — p _1 ||^|| 2 . With these definitions, equation (1) is 
equivalent to the following system for the vector (U, V) 



<u\ 




fu\_ 






{-A, 0) 







(39) 



The transformed wave field (u + , u ) and source (/ + , /_) were defined in (12) from 
(U, V) and (0, F). Recall that the matrix pseudodifferential operator Q is elliptic, 
and that Q^ 1 denotes a microlocal inverse, not an exact inverse, satisfying Q~ X Q ~ 

QQ- 1 rsj Id. 
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A technical complication is that the operators Q and Q^ 1 are not pseudodiffer- 
ential operators in (z, x, t) (only in (x, t)). We let ip = ip(z, x, D z , D x , D t ) be a 
microlocal cutoff around £ = r = 0. We let its symbol ip(z, x, C, £, t) be 1 for 
|C| < 2C|t|,|C| > 1 and for |C| > 3C|r|. Equation (13) is true if and only if 



d_ 

dz 



( 



p 1 



0. 



= ipQ 1 , microlocally on I 02 . (40) 



By [6, theorem 18.1.35] the operator ipQ _1 is a pseudodifferential operator with 
symbol that equals Q^ 1 modulo S~°° on Ig 2 . Using this theorem again we can see 
that in (40) a factor QQ~ l can be inserted in the left side before (U, V). So (40) 
holds if and only if 



ijQ- 1 



d_ 

dz 



p 



-A p / J 




microlocally on Ig 2 . (41) 



Therefore Theorem 1 follows from the lemma that we now state concerning the 
diagonalization of the operator 



Q- 1 



d_ 

dz 



( 



p 



-A p 0. 



Q. 



(42) 



Lemma 3 For suitably chosen Q and B, the operator (42) is equal to 



d_ 

dz 



B + 
B_ 



\ 



+ R, 



(43) 



where R = R(z, x, D x , D t ) is a 2 x 2 matrix pseudodifferential operator of order 1 
with symbol that is in S~°° on I' d2 . Here we can choose B± selfadjoint, with Q, B± 
satisfying (17) and (16). We can also choose Q with <5i,i = <5i,2 = 1, with B± 
satisfying (18) and Q satisfying 



Q 



i i \ 

-isgn(r)p^^a p ^ + order(O) i sgn(r)p~^a p ^ + order(O)^ 



(44) 



The proof mostly follows an argument of Taylor [8]. Some extra work is required 
to obtain the explicit expressions and the symmetry (self-adjointness) property. It 
is an order by order construction resulting in an asymptotic sum that is well defined 
according to standard results (see e.g. [6, Proposition 18.1.3]). The main tool is 
the composition formula, which says that the product of two pseudodifferential op- 
erators A(y, D y ) and B(y, D y ) is again a pseudodifferential operator with symbol 
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Aj^B given by the asymptotic sum 

T,-L- l ^A(y,r ] )d^B(y,r ] ). (45) 

a CX. 1 

Here a! = a^.-.aj (ify G R n ). 

For the computations in the proof of Lemma 3, we use pseudodifferential operators 
that are microlocally the square root or certain other powers of the operators A, 
A p . In addition we use an operator S(D t ) with symbol — sgn(r). In the following 
lemma we collect the needed information about these. 

Lemma 4 There is a pseudodifferential square root operator B = B(z, x, D x , D t ), 
satisfying B 2 ~ A microlocally on I' e . Its symbol satisfies 

n_1 db db 

B(z,x,£,t) = b + lib' 1 J2 J^JT + order (-l), on 4. (46) 

j=i °v c>x j 

There are fourth and second roots A p s , s = |, |, microlocally on I' d2 , with prin- 
cipal symbol that equals a p (z, x, £, t) s on I' e . There are inverse A p ~ s , s = |, |, 1 
microlocally on Ig 2 , with principal symbol a p (z, x, £, t)~ s on I' d2 . The operators 
B, A p s can be chosen selfadjoint. Let S(z, x, D x , D t ) be a pseudodifferential oper- 
ator with symbol equal to — sgn(r) on I' d2 , |r| > 1, and selfadjoint. Then B ~ SA2 
on I' e . We have 

B = Sp^A p ^p^ + order(-l) 

= Sp'^A p ^pA p ^p-^ + order(-l), on I' 0a . (47) 



PROOF. Following the standard argument (compare e.g. Lemma II. 6. 2 in [9]) we 
look for the square root as an asymptotic sum J2JLo with each a pseu- 
dodifferential operator of order 1 — j. We let have principal symbol b, then 

( T (o))2 = A + R (o) t with ^(0) of order 1 Now SU pp 0S e we have T®, j = 0,...,k, 

such that 

k \ 2 

2 T {j) =A + R {k) , with R {k) of order 1 - k. (48) 
3=0 J 

Then we let T^ k+1 ^ have principal symbol — ^b~ 1 R^ k \ and (48) is valid with k + 1 
instead of k (note that this choice of T^ k+1 ^> is unique to highest order on I' e2 ). In 
this case is i YTjZl |r J| + order(-l), which leads to (46). Because b is real 

all the can be chosen selfadjoint. The square root A} follows similarly and its 
fourth root equals the square root of the square root. The existence of microlocal 
inverses is standard (see Theorem 18.1.9 in [6]). 

The symbol of S is locally constant on I' d2 so S commutes microlocally with p s and 
the A p . Also the commutator [A p s , p s '] is a lower order operator (of order | — 1). 
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To compute the square of Sp * A p 2 p 4 we commute a factor p 4 to the left of the first 
A p 2 and one to the right of the second A p z . Using that the multiple commutator is 
again an order lower, we find 

(Sp*A p *p*) 2 = p^App^ + order(O) = A + order(O), on I' 92 . (49) 

Equation (47) is clearly valid for the principal symbols. It follows from (49) and 
the fact that T^> is unique to highest order that the first equality in (47) is satisfied. 
The second equality of (47) follows similarly. □ 



Proof of Lemma 3 Commuting Q 1 and ^ we find that (42) is equal to 



d_ 

dz 



-Q- 1 



V 



>-^o- 

-A„0 dz 



(50) 



We first consider the second term of (50), which is a contribution of order 1 (the 
third term is of order 0). The eigenvalues and eigenvectors of the principal symbol 



matrix 



p 



-a p 0, 



are given by 



eigenvalues : ± ip(z, x) 1 ^ 2 a p (z, x, £, r) i/2 , eigenvectors : 



1/2 



We first make a highest order choice Q = such that the matrix Q 1 
becomes diagional on I' 02 . Below we will add lower order terms. We set 



p(z,xy/* 

(51! 

p 



-A p 0, 



Q 



Q 



(o) 



A -A 1 

A p 4 pi 



1 1 

)4 



iSApip 4 -iSAptp 4 



(52) 



with S as in Lemma 4. Outside (9) we still require that the symbol is an 
invertible matrix symbol of order 2 2 ] . It is easily seen that this is possible. 

\ 2 2 

The inverse of Q(°) satisfies 

Q (0) (^^^) _1 = 



/ - 1 A A . „ 1 , _I ' 
' p 4^ p 4 _ 1 5 , p 4 A j0 4 

1 . I 



P 

, _ 1 



\p ±Ap± \Sp*A p 4, 



, on I' e 



e 2 - 



(53) 
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With this choice we find, using the equalities (47) 



Q 1 I ° P \Q 



-A p 0, 




order(-l), on I' e 



(54) 



where B is the square root operator defined in (46). 

Next we take the third term in (50). It can be seen easily that 

dA 



3 A s 

-^A-° = s^A,- 1 + order(-l), on I', y 



(55) 



It follows that 



dQ~ l 

dz 



I 



Q 



\( d -^Ap - i -sep- i V 

4 \ dz P dz" ) 



+ order(-l), on I' d2 . 

(56) 



d_ 

dz 



B 



+ on/;,. (57) 



vKtSt^" 1 - ft"" 1 ) 

Thus, with Q = Q 1 - -*, the expression (50) is equal to 
/ 

Here R is a pseudodifferential operator of order 1, that is of order —1 on the set I' d2 . 

In expression (57) the highest, first order part is diagonal, while there are lower 
order off-diagonal terms. Following [8] we will remove the off-diagonal terms order 
by order. To remove zeroeth order off-diagonal terms we modify Q, and set it equal 
to Q = QM = Q^(l + KMy Here is an operator that remains to be chosen, 



is of order —1 and of the form KM 



( o <)' 

K2° 



This results in an additional 



contribution to (57) of order zero, which is given to highest, zeroeth order by 



\B 

-LB 



KM -KM 




\(BK$+K$B) 



\-\{BKM+K^B) 



(i) 







/ 

(58) 

It follows that there is a symbol RM such that this contribution cancels to ze- 
roeth order off-diagonal contribution on I' e , not changing the first and zeroeth 
order diagonal part. By considering further modifications of the form Qk> +1 ) = 
QM{\ + K^ +1 ^), of order —j, also the lower order off-diagonal terms can 

be removed. This proves the existence of B± satisfying (16). 
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Next we prove the self adjointness. The operator satisfies 



1 \ „ / -iS 

Q (0) * 

0-1 us 



Q 



(0)-l 



(59) 



It follows that with Q = we have 



Q- 1 



o P 



-A p 0, 



Q = - 




Q' 1 



o p 1 



g 




(60) 



Hence with Q = the second term (50) is the sum of an anti-selfadjoint diag- 
onal part and a selfadjoint off-diagonal part (this also follows from explicit com- 
putation). The same property is true for 9 ® z Q. To prove that B± can be chosen 
selfadjoint it is sufficient to show that this property is still true when the Q is mod- 
ified order by order as described above. So suppose the property holds for 
Then there is a off-diagonal, self-adjoint K^ +1 \ with correct highest order term 
as above. We modify the definition of Qk> +1 ) with lower order terms according to 
QU+ 1 ) = qU) exp(K^ +1 ^) (meaning the power series for exp). Then the microlo- 
cal inverse is given by Q&+ 1 )- 1 = exp(— K^ +l ^)Q^^ 1 . It is easy to see from the 
power series for exp that 



(x o^ 



exp(-i^' +1) ) 



1 
-1 



exp(i^ +1) ) 



(61) 



It follows that the new matrix valued differential operator stays the sum of an anti- 
selfadjoint diagonal part and a selfadjoint off-diagonal part. This shows the selfad- 
jointness. 

To compute the second choice of Q and B± we replace Q with QD where 



£) = diag(Qli 1 ,Qr 2 1 ) 



(62) 



If we denote the second choice of Q by Q, it follows easily that 



Q- 1 



d_ 

dz 



V 



p 1 



AO, 



d . I B + ,9,0- 

Q = Tz- 1 \ o B _ ^ 



-D + R. 



(63) 



The operator D is diagonal and the principal symbol of D is given by diag(p 2d*, 
p _ 2 a*). This shows the last part of the lemma. □ 
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4 Proof of Theorem 2 



We prove only the + sign, the — sign proceeds in the same way. Define v + , t>_ by 

(:)-'(") 

Because of (23) it follows that V-(z ) = on I' 02 fl {z = z }, hence 

v + (z ) = Q^U(z ). (65) 

Since on J+(z , 9 2 ) we have — r -1 ^ > 0, it follows that U = Q+v + on J+(z 0: 9 2 ). 
Equation (22) implies that 

for some / satisfying / = on Ig 2 . Let w+ be the solution to 

P ,+w = 0, w+\ Z0 = v + \ ZQ . (67) 

By the initial condition and a propagation of singularities result, and the fact that 
P+(v + — w + ) ee on Iq 2 , it follows that 

w+ ee w + + 5-, (68) 

where g = Oon J + (z ,^). Since the symbol of is in S~°° outside J + (z , 9 2 ), 
we have K + g = 0. Hence 

U + EE EE X + K - (?) EE K + (Q^U). (69) 

Therefore 

g+M+EE^ + tg+^jg; 1 )^. (70) 

For the commutator term [Q + , K]Q^ it follows from (34) and (35) that 

[Q^KjQ- 1 ~Oon + J(z o ,0 1 ), (71) 
[Q+, K]Q- 1 e Op outside J + (z , 9 2 ), if z - z > 5. (72) 

This completes the proof of the theorem. 
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